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ABSTRACT 

Magnetic fields of low-mass stars and planets are thought to originate from 
self-excited dynamo action in their convective interiors. Observations reveal a 
variety of field topologies ranging from large-scale, axial dipole to more structured 
magnetic fields. In this article, we investigate more than 70 three-dimensional, 
self-consistent dynamo models obtained by direct numerical simulations. The 
control parameters, the aspect ratio and the mechanical boundary conditions 
have been varied to build up this sample of models. Both, strongly dipolar and 
multipolar models have been obtained. We show that these dynamo regimes can 
in general be distinguished by the ratio of a typical convective length-scale to 
the Rossby radius. Models with a predominantly dipolar magnetic field were 
obtained, if the convective length scale is at least an order of magnitude larger 
than the Rossby radius. Moreover, we highlight the role of the strong shear 
associated with the geostrophic zonal flow for models with stress-free boundary 
conditions. In this case the above transition disappears and is replaced by a region 
of bistability for which dipolar and multipolar dynamos co-exist. We interpret 
our results in terms of dynamo eigenmodes using the so-called test field method. 
We can thus show that models in the dipolar regime are characterized by an 
isolated 'single mode'. Competing overtones become significant as the boundary 
to multipolar dynamos is approached. We discuss how these findings relate to 
previous models and to observations. 
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Introduction 



Topologies of stellar and planetary magnetic fields stemming from hydrodynamic dy- 
namo action are highly variable. Observations revealed complicated field configurations 
dominated by higher multipol es, but also large-scale or even dipolar stellar and planetary 



magnetic fields are ubiquitous ([Anderson et al.ll2010l : iHulot et al. 



2OIOI : iRussell fc Doughertv 



2010 ; Donati fc Landstreet 20091 ). The latter is the more astonishing, as the magnetic 



Reynolds number in planetary and stellar convection zones is often large and most hydrody- 
namic dynamos operate far above their threshold. Observations and numerical simulations 
indicate that rapid global rotation and thus the ordering influence of the Coriolis force is 



2004; 


Kapvla et al. 


2009; 


Brown et al. 


2010 


)■ 



the ex istence of a dipolar and a multipolar dynamo regime and IChristensen fc Aubert 



( 120061 ) showed that the transition between both is governed by a local Rossby number, 
i.e. the influence of inertia relative to the Coriolis force. Similar results were reported by 



Sreenivasan &: JonesI (120061 ). too. Dipolar models were found for small Rossby numbers; they 
are separated by a fairly sharp regime boundary from multipolar models, for which inertia 
is more important. 



Sreenivasan fc JonesI (120111 ) argued that dipolar magnetic fields enhance the kinetic he- 
licity and are therefor e easier to maintain than fie lds with a more complicated field topology. 
However, as noted by ISreenivasan fc JonesI (1201 ll ). the r elation between kineti c helicity and 
induction mechanisms is not straightforward. Moreover, ISchrinner et al.l ( 120071 ) showed that 
the kinetic helicity is indeed a bad proxy measure for the induction effects (a-effect) in these 
models. 



Schrinner et al.l ( l2010a( ) pointed out that the high relative dipole field strength for mod- 
els in the low Rossby number regime is associated with the dominance of only one dipolar 
eigenmode. Spatially more complex modes were found to be strongly diffusive. In this 
dynamo regime, higher order contributions to the magnetic field result mainly from the de- 
formation of the fundamental mode by th e turbulent flow. Subsequently, th e small scale 
contributions decay due to ohmic diffusion (lHoyngll2009l ; ISchrinner et al.l 1201 ll ). As a result 
of the dominance of only one, real eigenmode, the axis of the dipole field is stable, polarity 
reversals or oscillations of the magnetic field do not occur. 

However, the reason for the dominance of an isolated eigenmode at low Rossby numbers 
and the cause of the dipole breakdown with decreasing infiuence of the Coriolis force are 
at present not well understood. In this study, we investigate how rotation infiuences the 
field topology in dynamo models in some more detail. We show that rotation-dominated 
convection in a spherical shell leads to a distinctive azimuthal field morphology which is well 
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represented by the fundamental dynamo mode, but very different from the field pattern of 
the following higher-order modes. We argue that this discrepancy is responsible for the clear 
preference of the fundamental mode over higher-order modes. Furthermore, we discuss how 
the dominance of the fundamental mode is reduced if convection becomes less affected by 
rotation and the typical length scale of convection decreases relative to the Rossby radius. 
Because the topology and the time variability of the magnetic field in numerical dynamo 
simulations are closely correlated we likewise address the question under which conditions 
models exhibit fairly coherent oscillations, irregular polarity reversals or a stable dipole field. 

Our analysis is based on 72 dynamo models in the Boussinesq approximation with dif- 
ferent aspect ratios and me chanical boundary con d itions . We generalize here the Rossby 
number criterion given by IChristensen fc AubertI (120061 ) to our sample of models. For 
some of them, we com pute mean-field coefficients with the help of the test-field method 
( ISchrinner et al.l 120071 ) to reveal their dynamo mechanisms. 



2. Dynamo Calculations 

Our dynamo models are solutions of the MHD-equations for a conducting Boussinesq 
fluid in a rotating spherical shell. The fluid motion is driven by convection due to an 
imposed temperature difference, AT, between the inner and the outer shell boundaries. The 
fundamental length scale of our models is the shell width L, we scale time b y 1? jv, with v the 



kinematic viscosity, and temperature is scaled by AT. Moreover, following IChristensen et al. 



( 120011 ). the magnetic field is considered in units of y/gjrrjn, with g denoting the density, 
the magnetic permeability, t] the magnetic diffusivity and Q the rotation rate. With these 
units, the dimensionless momentum, temperature and induction equation are 



E {— + v - Vv - V^v j + 2z X V + VP = 

Ra—T + -^{VxB)xB, (1) 

To Fm 

dT 1 , 

- + «.VT=-V=r. (2) 

_ = Vx(.xB) + — V^B. (3) 

We also note that v and B are solenoidal. The system of equations is governed by four 
dimensionless parameters, the Ekman number E = u/QL"^, the (modified) Rayleigh number 
Ra = axgo^TL/uQ, the Prandtl number Pr = zz/k, and the magnetic Prandtl number 
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Pm = v/rj. In these definitions, ut stands for tlie tliermal expansion coefficient, go is tlie 
gravitational acceleration at the outer boundary, and k is the thermal diffusivity. A further 
control parameter is the aspect ratio of the shell defined as the ratio of the inner to the outer 
shell radius, x = ^i/^o- It determines the width of the convection zone. 

The mechanical boundary conditions are either i) no slip at the inner and the outer 
boundary, ii) no slip at the inner and stress free at the outer boundary or iii) stress free at 
both boundaries. In the latter case, the angular momentum in the direction of the rotation 
axis was conserved. Furthermore, the magnetic field matches a potential field outside the 
fluid shell and fixed temperatures are prescribed at both boundaries. 

Some of the models investigated here exhibit bistability, i.e. the solution realized de- 
pends on the initial conditions for the magnetic field. We started most of the simulations with 
a dipolar magnetic field, but varied its initial amplitude over several orders of magnitude. 
Some calculations were started from another model with slightly different parameters to test 
for hysteresis. In the bistable regime, models resulting from simulations with an initially 
weak magnetic field are here referred to as the 'non-dipolar branch'. They are distinguished 
from dipolar solutions initially started with a strong magnetic field. 



The numerical solver used to compute solutions of equations ([I])-© is PaRoDy (jPormy et al 



and furthe r developments). The numerical method is similar to that described in 



Glatzmaierl (|1984( ) except for the radial discretisation, which is treated in physical space on 
a stretched grid (allowing for a parallelization by a radial domain decomposition). More- 
over, the pressure term has been eliminated by considering the double curl of the momentum 
equation. 



3. Non-Dimensional Output Parameters 



Our numerical dynamo-models are characterized by non-dimensional output parameters. 
Dimensionless measures for the flow velocity are the magnetic Reynolds number, Rm = 
VT^tiis L/r], and the Rossby number, Ro = Vrms/^L. In both definitions, frms stands for the 
rms velocity of the flow. Similarly, -Brms denotes the rms field strength. We also introduce a 
local Rossby number, Ro£ = Ro ip/rf, based on the mean harmonic degree ip of the poloidal 
velocity field, 

,< Me ■ Me > 



E 



<Vp-Vp> 



(4) 



The brac kets in (HI) denote an average o ver time and radii. In contrast to the definition intro- 
duced by IChristensen fc AubertI (120061 ). our modified Rossby number relies on a convective 
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length-scale and not on the mean half wavelength of the total flow (see also the Appendix). 

The magnetic field strength is measured by the dimensionless Lorentz number, Lo = 
Brms/ (y/oJ^^L), and the classical Elsasse r number A = B'^^^/flp u p. Th ey are related through 



A = Lo^Pm/E. Moreover, following IChristensen fc AubertI ( l2006l ). we characterize the 



geometry of the magnetic field by the relative dipole field strength, /dip, which is defined as 
the ratio of the average field strength of the dipole field to the field strength in harmonic 
degrees £ = !,■■■ , 12 at the outer boundary. 

A non-dimensional measure for the heat transport is given by the Nusselt number, Nu, 
defined as the ratio of the total heat flow to the conducted heat flow, Qcond = 4:TTrarigcKAT/ L 
with the heat capacity c. 



4. Mean-Field Analysis 

We present dynamo mode ls obtained by direct riumerical simulations (DNS) and do not 
use the mean- field formalism ( ISteenbeck et al.l Il966l : iMoffattI Il978t iKrause &: Radlerl Il980[ ) 
to predict any dynamo action. However, the mean-field approach provides useful theoretical 
concepts and mathematical tools to interpret numerical dynamo models. It is usually set up 
by splitting the velocity and the magnetic field in mean and fluctuating parts, v = V + v' 
and B = B + b'. Mean quantities are here denoted by an overbar and defined as azimuthal 
averages. The objective of mean-field theory is to predict the evolution of the averaged or 
mean magnetic field, B, which is in general determined by the dynamo equation. 



dB 



V X {£ + V X B - riV X B). 



(5) 



The crucial new term in (E]) compared to the induction equation is the mean electromotive 
force, £ = v' X h' . It is a functional of d, V, and B, which is affin-linear in B. Provided that 
there is no small-scale dynamo action, £ is also homogeneous in B. If, moreover, the mean 
electromotive force depends only instantaneously and locally on B, it may be parameterized 
in terms of the mean field and its first derivatives, 



cxB -J X B - /3(V X B)-d X (y X B)- kVB. 



(6) 



The parameters, also known as mean-field coefficients, are vectors (7 and d), symmetric 
tensors of second rank [ex. and f3) and a third rank tensor (k). They depend only on the 
velocity field and the magnetic diffusivity, but not (explicitl y) on the mag netic field. A 
physical interpretation of the mean- field c oefficients is g iven by iRadlerl (119951 ). The a-tensor 
describes the classical a— effect (see also iParkerl Il955l ). 7 is associated with an advective 
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transport of the mean magnetic field, /3 can be interpreted as a turb ulent diffusivity , 5 
may contribute to a shear-current effect, an inductive effect first noted by iRadler f ll969al lb[). 
whereas the K-terms are more difficult to interpret. The altogether 27 independent mean- 
field coefficients w ere determined for several ri umerical models with the help of the so-called 
test-field method JSchrinner et al.lboosl . [2007I ). 



A comparison between DNS and mean-field calculation s revealed that parameterization 
(ISl) is indeed reliable for a wide class of dynamo models (jSchrinneii l2011t ISchrinner et al. 
201X1 ). An important limitation is the intrinsic kinematic character of the mean-field ap- 
proach followed here. In general, it is therefore only applicable to dynamo models which 
can be repr oduced kinemati c ally an d belong to the so-called kinematically stable regime 
identified by ISchrinner et al.l ( l2010a( ) . 



For a mean-field analysis of some of our dynamo models derived from ([I])-®, we solve 
the dynamo equation as an eigenvalue problem 

aB = V X DB , (7) 

with the time-averaged dynamo operator D defined through 

DB = vxB- aB -jxB- f3{V xB)-Sx{VxB)- kVB. (8) 

On average, the evolution of the mean field is then proportional to exp (o"jt) with eigen- 
vectors and e igenval ues CTj. More details about the eigenvalue calculation are given in 
Schrinner et al.l ( l2010bl ). 



5. Results 

We present results from 72 dynamo models obtained by varying all four dimensionless 
control parameters, the aspect ratio of the shell as well as the mechanical boundary condi- 
tions. Table [1] is divided in three sections for three different types of mechanical boundary 
conditions and provides the control parameters and the output parameters introduced above 
for each model. Within each section, the models are listed in the order of increasing lo- 
cal Rossby number. For some models with stress-free mechanical boundary conditions, two 
different solutions coexist depending on the initial conditions for the magnetic field. Thus, 
these models appear in pairs and are labeled by the letter 'd' (dipolar) if the initial magnetic 
field was strong and by 'm' (multipolar) if the simulations were started from a weak magnetic 
field. 



Kutzner fc ChristensenI (120021 ) identified a dipolar and a multipolar dynamo regime. 
Both regimes can be recovered here if the models are ordered by the local Rossby number 
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introduced in section |3] (see also IChristensen fc Aubertll2006l ). Figure [T] (left panel) shows 



two distinct dynamo regimes: models with a dipole dominated magnetic field at low Rossby 
numbers (filled symbols) and models with a more multipolar magnetic field at higher Rossby 
numbers (open symbols). A fairly sharp transition between both regimes occurs at Roi ~ 0.1. 
The exact value may depend weakly on the choice of the mechanical boundary conditions and 
seems to be closer to Roe ~ 0.12 for models with rigid boundaries. The right panel of Fig. 
[1] illustrates that models of both regimes differ in their field topology, but not necessarily 
in their field strength. We emphasize that only the definition of a local Rossby number 
based on a convective length scale (i.e. not taking into account the zonal fiow) enables us to 
extend the well known dipole-multipole partition to models with different aspect ratios and 
boundary conditions. This is most clearly demonstrated by the sequence of models 29, 31, 
32, 34 and 35 from Tabled] which is depicted in Fig. [21 For these models, the Ekman number, 
the Rayleigh number (normalized by its critical value), and the Prandtl numbers were kept 
constant whereas the aspect ratio was progressively increased. The thinner the convection 
zone, the smaller the convective length scales, and the mean harmonic degree of the poloidal 
velocity field increases systematically from ip = 11 for model 29 with x = 0.5 to = 20 
for model 35 with x = 0.65. Consequently, Roe grows from Roe = 4.70 x 10^^ for model 29 
to Roi = 1.14 X 10~^ for model 35, which already belongs to the multipolar regime. The 
decrease of the convective length scales caused a transition towards the multipolar regime 
and is adequately measured by the local Rossby number introduced here. If, however, Ip 
were derived from the total velocity field, neither ip nor Rog would increase in the same 
way, owing to a major contribution of the mean zonal fiow to the kinetic energy density in 
model 35. The Rossby number criterion would then fail to predict the observed transition 
between both dynamo regimes. Note that decreasing the Ekman number from E = 10~^ 
to E = 3 10""^ in model 30 leads again to a dipolar field despite the high aspect ratio of 
X = 0.65 and that we find a multipolar magnetic field for model 37 at lower aspect ratio, 
X = 0.6, but higher Rayleigh number. Both examples illustrate that the transition is indeed 
characterized by Roe. However, the definition of Roe is empirically motivated, there may be 
other ways to define a relevant Rossby number to distinguish both dynamo regimes; some 
of them are discussed in Appendix A. 

Moreover, Fig. [3] demonstrates that the relative dipole field strength does not simply 
depend on the magnetic Reynolds number. The distance from the dynamo threshold of a 
model does not determine its field topology. 

An exception to the Rossby number criterion for the dipolarity of the magnetic field in 
dynamo models occurs in the case of stress free boundary conditions. Simulations started 
from a weak magnetic field (m-models) are not simply dipole dominated, independent of 
their local Rossby number. These models left out in Fig. [1] were included in Fig. H] and 
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Fig. 1. — Left: Relative dipole field strength versus the local Rossby number for all models 
from Table [T] apart from those with stress-free boundary conditions belonging to the multi- 
polar branch (m- models). Filled symbols stand for models dominated by a dipole field, open 
symbols denote multipolar models. The symbol shape indicates different types of mechan- 
ical boundary conditions: circles mean no-slip conditions at both boundaries, triangles are 
models with a rigid inner and and a stress-free outer boundary, and squares stand for models 
with stress-free conditions at both boundaries. Right: Elsasser number versus local Rossby 
number for the same sample of models. 
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are represented by diamonds. As apparent from Fig. |H dynamo models with stress free 
mechanical boundary conditions form two branches, an upper branch of dipolar models for 
Rop < 0.1 and a lowe r branch of models wit h a less constrained field geometry (see also 
Busse fc SimitevI l2009t ISimitev fc Bussd |2009[ ) . We refer to them as the 'dipolar' and the 
'multipolar' branch instead of the 'weak-field' and the 'strong-field' branch to avoid confusion 
with corresponding notions introduced in the limit of vanishing viscosity and inertia in 
the context of planetary dynamos. Nevertheless, the saturated magnetic field strength, as 
measured by the Lorentz number, is always larger for dipolar than for multipolar models. 
For Rog > 0.1, both branches coincide and their distinction is no longer meaningful (see 
Fig JT^ . The region of bistability in parameter space, however, do es not solely depend o n 
the local Rossby number, but also on the magnetic Prandtl number (IBusse fc Simitevll2009l ): 
for models 40, 42, and 44, a weak-field branch does not exist. 

There is a strong correlation between the topology and the time dependence of the 
magnetic field in dynamo models. Sudden polarity reversals or oscillations of the magnetic 
field do not occur in dipole dominated models in the low Rossby number regime. Conversely, 
reversals and oscillations are frequent in non-dipolar models with Ro£ > 0.1 as well as in 
models with lower local Rossby numbers with stress-free boundary conditions belonging to 
the multipolar branch. Whether non-dipolar models exhibit fairly coherent oscillations or 
irregular reversals of the magnetic field strongly depends on the magnetic Reynolds number. 
Coherent oscillatory solutions of the induction equation are most clearly visible in so-called 
butterfly diagrams; contours of the azimuthally averaged radial magnetic field at the outer 
boundary are plotted versus time and colatitude. Figure O gives two examples. The left 
panel shows a very coherent dynamo wave at Rm = 102 (model 45m), whereas the but- 
terfly diagram at Rm = 258 on the right-hand side (model 38) is much less periodic and 
a cycle period cannot be identified. Dynamo models (in the non-dipolar regime) at higher 
magnetic Reynolds number exhibit even less temporal coherence. Following this somewhat 
arbitrary and qualitative criterion, we find that non-dipolar dynamos of our sample with 
Rm < 200 generate magnetic fields which vary periodically in time. The lower the magnetic 
Reynolds number, the more coherent is the time variability of the magnetic field. A third 
property is intimately l inked to the topology and the time dependence of the magnetic field. 
Schrinner et al.l ( l2010a( ) looked at the evolution of a passive vector field described by a second 
induction equation which was solved simultaneously with equations ([I])-(l3]). The magnetic 
'tracer-field' was advanced at each time-step with the self-consistently determined velocity 
field but did not con t ribut e to t he Lorentz force. A simil a r expe rimen t had been perfo r med by 
Cattaneo fc TobiasI ( l2009l ) and lTilgner fc Brandenburg ( l2008l ). too. ISchrinner et al.l ( l2010al ) 
found that the tracer field grows exponentially for multipolar and reversing models, whereas 
it remains kinematically stable for dipolar dynamos in the low Rossby number regime. This 
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Fig. 2. — The relative dipole field strength for models 29, 31, 32, 34 (filled triangles) and 
model 35 (open triangle) and the mean harmonic degree ip of the poloidal velocity field (right 
panel) versus the local Rossby number. The aspect ratio (upper x-axis) has been gradually 
increased for this sequence of models leading to larger ip and thus also to larger Roe. Model 
35 has undergone a transition from the dipolar to the multipolar dynamo regime and the 
dipole field strength has dropped drastically. 
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Fig. 3. — The magnetic Reynolds number versus the local Rossby number for all models 
from Table [TJ For the meaning of the symbol style we refer to the caption of Fig. [TJ 




Fig. 4. — The relative dipole field strength versus the local Rossby number as in Fig. [T]but 
here with the m- models included (diamonds). All other notations as in Fig. [T] 



relation can be extended to the sample of models studied here and serves to distinguish both 
dynamo classes. Dipolar models represented by filled symbols in Fig. H] are kinematically 
stable, all multipolar models represented by open symbols are kinematically unstable. 
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Table 1. Overview of the simulations carried out, ordered with respect to their modified 
Rossby number and their mechanical boundary conditions. 
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89 


6.63 X 10-3 


36 


1 X 10-3 


150 


5 


1 


0.60 


16 


1.17 X 10-1 


0.83 


113 


4.09 X 10-2 


37 


1 X 10-3 


175 


5 


1 


0.60 


16 


2.08 X 10-1 


0.18 


205 


1.71 X 10-2 


38 


1 X 10-3 


200 


5 


1 


0.60 


16 


2.68 X 10-1 


0.20 


258 


1.87 X 10-2 


39 


1 X 10-3 


250 


5 


1 


0.60 


16 


3.40 X 10-1 


0.12 


331 


2.55 X 10-2 








Stress- Free Boundary 


Conditions 








40 


1 X lO"* 


365 


2 


1 


0.,35 


12 


1.82 X 10-2 


0.92 


92 


5.69 X 10-3 


41d 


3 X 10-5 


600 


1 


1 


0.35 


16 


2.11 X 10-2 


0.86 


140 


3.82 X 10-3 


41m 


3 X 10-5 


600 


1 


1 


0.35 


16 


2.56 X 10-2 


0.26 


144 


2.35 X 10-3 


42 


1 X 10-* 


375 


4 


1 


0.35 


13 


2.20 X 10-2 


0.70 


216 


8.80 X 10-3 
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Table 1 — Continued 



Model 






E 


Ra 


Pm 


Pr 


X 




Ro 


I 




/dip 


Rm 


Lo 


43d 


1 


X 


10-4 


462 


2 


1 


0.35 


16 


3.04 


X 


10- 


. 2 


0.84 


121 


8.44 


X 10-3 


43m 


1 


X 


10-4 


462 


2 


1 


0.35 


16 


3.49 


X 


10- 


. 2 


0.33 


146 


4.91 


X 10-3 


44 


1 


X 


10-4 


750 


6 


1 


0.35 


17 


4.20 


X 


10" 


. 2 


0.77 


462 


2.50 


X 10-2 


45d 


1 


X 


10-4 


586 


1 


1 


0.35 


18 


4.54 


X 


10" 


. 2 


0.93 


78 


9.71 


X 10-3 


45m 


1 


X 


10-4 


586 


1 


1 


0.35 


18 


5.73 


X 


10" 


. 2 


0.20 


102 


5.60 


X 10-3 


46d 


1 


X 


10-4 


749 


2 


1 


0.35 


18 


5.30 


X 


10" 


. 2 


0.88 


187 


1.40 


X 10-2 


46m 


1 


X 


10-4 


749 


2 


1 


0.35 


18 


7.00 


X 


10" 


- 2 


0.39 


244 


8.12 


X 10-3 


47d 


1 


X 


10-4 


750 


4 


1 


0.35 


17 


4.20 


X 


10" 


■2 


0.85 


320 


2.06 


X 10-2 


47m 


1 


X 


10-4 


750 


4 


1 


0.35 


20 


7.20 


X 


10" 


-2 


0.29 


460 


9.89 


X 10-3 


48d 


1 


X 


10-4 


750 


0.75 


1 


0.35 


19 


6.03 


X 


10" 


■2 


0.95 


77 


1.11 


X 10-2 


48m 


1 


X 


10-4 


750 


0.75 


1 


0.35 


18 


7.70 


X 


10" 


■2 


0.20 


100 


7.02 


X 10-3 


49d 


3 


X 


10-4 


510 


2 


1 


0.35 


13 


7.26 


X 


10" 


-2 


0.93 


120 


2.08 


X 10-2 


49m 


3 


X 


10-4 


510 


2 


1 


0.35 


12 


9.15 


X 


10" 


-2 


0.25 


155 


1.09 


X 10-2 


50d 


1 


X 


10-4 


750 


0.5 


1 


0.35 


12 


7.60 


X 


10" 


-2 


0.97 


53 


1.03 


X 10-2 


50m 


1 


X 


10-4 


750 


0.5 


1 


0.35 


12 


8.20 


X 


10" 


-2 


0.21 


70 


6.32 


X 10-3 


51d 


1 


X 


10-4 


750 


2 


0.3 


0.35 


15 


7.90 


X 


10- 


-2 


0.86 


340 


2.55 


X 10-2 


51m 


1 


X 


10-4 


750 


2 


0.3 


0.35 


15 


1.06 


X 


10" 


1 


0.38 


450 


1.52 


X 10-2 


52d 


1 


X 


10-4 


1110 


1 


1 


0.35 


20 


9.20 


X 


10- 


-2 


0.92 


146 


1.61 


X 10-2 


52m 


1 


X 


10-4 


1110 


1 


1 


0.35 


20 


1.19 


X 


10" 


1 


0.33 


188 


1.13 


X 10-2 


53 


3 


X 


10-4 


510 


1.5 


1 


0.35 


12 


9.46 


X 


10" 


-2 


0.20 


122 


1.01 


X 10-2 


54 


1 


X 


10-4 


1333 


2 


0.3 


0.35 


15 


1.60 


X 


10" 


1 


0.38 


450 


1.52 


X 10-2 


55 


1 


X 


10-4 


1000 


1 


1 


0.35 


20 


1.13 


X 


10" 


1 


0.16 


176 


9.83 


X 10-3 


56 


1 


X 


10-4 


1200 


1 


1 


0.35 


21 


1.36 


X 


10" 


1 


0.16 


204 


1.18 


X 10-2 


57 


1 


X 


10-4 


1372 


1 


1 


0.35 


19 


1.38 


X 


10" 


1 


0.30 


234 


1.32 


X 10-2 


58 


1 


X 


10-4 


1757 


1 


0.3 


0.35 


15 


1.60 


X 


10" 


-1 


0.38 


330 


1.61 


X 10-2 


59 


1 


X 


10-4 


1497 


1 


1 


0.35 


21 


1.67 


X 


10" 


-1 


0.19 


245 


1.45 


X 10-2 


60 


1 


X 


10-4 


1627 


1 


1 


0.35 


22 


1.82 


X 


10" 


-1 


0.23 


266 


1.52 


X 10-2 


61 


1 


X 


10-4 


1823 


0.3 


1 


0.35 


22 


2.35 


X 


10" 


-1 


0.19 


101 


1.45 


X 10-2 


62 


3 


X 


10-4 


2000 


2 


0.3 


0.35 


15 


4.81 


X 


10" 


-1 


0.15 


658 


3.63 


X 10-2 



Note. — Some of the models in Table [T] were already published in ISchrinner et al.l l l2010al ') with a 
somewhat different definition of £p and thus also of Roi. 
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0.1 0,2 0,15 0,3 

Magnetic diflueion time Magnetic diflueion time 

Fig. 5. — Left panel: Butterfly diagram for model 45m, Rm = 102. Right panel: Butterfly 
diagram for model 38, Rm = 258. The contour plots are normalized by their maximum 
absolute value and the grey-scale coding varies from +1 (black) to -1 (white). 



E 




0.01 



0.10 

Ro, 



1.00 



Fig. 6. — Magnetic Reynolds number versus the local Rossby number for all multipolar 
models of Table [H A dot in an open symbol indicates that a coherent dynamo wave was 
found. The meaning of all other notations is as defined in Fig. [1] and Fig. HI 
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6. Discussion 

6.1. Reason for the dipole breakdown 

What causes the dipolarity of dynamo models far above the dynamo threshold in the 
low Rossby number regime? The Rossby number criterion supported by this work predicts a 
dipole dominated magnetic field, if the typical length scale of convection, (^p/vr) L, is at least 
one order of magnitude larger than the Rossby radius, i.e. 0.1 (?p/7r) L > Wrms/^- Convection 
is then strongly infiuenced by rotation and organized in columns parallel to the rotation axis. 
A mean-field analysis reveals that the magnetic field generated and maintained by these con- 
vective rolls is dominated by only one real dipolar eigenmode with approximately zero growth 
rate. All more structured and in general complex overtones (i.e. higher order modes) are 
highly diffusive. The kinematic stability of models at low Roi as we ll as their comparativel y 



steady dipole field are a consequence of this single mode property (jSchrinner et al.ll2010af ). 
However, the reason for the breakdown of the dipole field at Roe ^ 0.1 is at present not well 
understood. We computed the eigenvalue spectra and the eigenmodes of the time-averaged 
dynamo operator D for a sequence of models with increasing Ro£ to gain more insight. The 
sequence consists of model 29, models 31-32, and model 34, already introduced in section |5l 
These are kinematically stable models {Roe < 0.1) with mixed mechanical boundary condi- 
tions and aspect ratios varying from x = 0.5 to x = 0.63. Models with the same parameter 
values and boundary conditions but an aspect ratio lower than 0.5 do not exhibit any dy- 
namo action. We therefore considered in addition model 4 with rigid boundary conditions 
as an example of a dynamo model at lower Rossby number. 

Figure [7] shows the growth rates of the fundamental mode and of the first two dipolar 
overtones versus the local Rossby number for these five models. The fundamental modes have 
on time average approximately zero growth rate, as it is expected for saturated dynamos, 
all overtones are diffusive. While there is typically a large gap between the fundamental 
mode and the first overtone for models at low Roe, both growth rates are much closer if 
the Rossby number increases. Quadrupolar modes were omitted in Fig. [7| because the field 
realized in the DNS is of purely dipolar symmetry for the five examples considered. For 
completeness, the eigenvalues of the first dipolar and quadrupolar modes are listed in Table 
[2J We emphasize again that the change in the spectra visible in Fig. [7] is not correlated to 
an increase of Rm. The magnetic Reynolds number is highest for model 4 [Rm = 86) and 
does not change significantly for the following models at larger Roe. 

At low Rossby numbers, the fundamental mode is well separated from the following 
modes (or overtones) which are strongly damped (see Fig. [7]). We refer to this characteristic 
as the 'single-mode' property. As the Rossby number is increased, the average growth rates 
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Fig. 7. — Growth rates of the fundamental mode (triangle) and of the first two dipolar 
overtones (stars) for five models at different Rossby number; in the order of increasing 
Rossby number, these models are: model 4, model 29, model 31, model 32, and model 34. 



Table 2. Eigenvalues for the first modes with dipolar (a") and quadrupolar (a^) 
symmetry for the five models considered in Fig. [7] in units of rj/L'^. 



Model 






<^2 


'^0 




4 


(-3.87,0.00) 


(-34.8, ±10.3) 


(-40.6,0.00) 


(-6.30,0.00) 


(-28.7, ±2.69) 


29 


(0.50,0.00) 


(-16.1,0.00) 


(-25.5, ±4.21) 


(-18.3,0.00) 


(-25.5, ±6.19) 


31 


(3.30,0.00) 


(-11.1,0.00) 


(-18.8, ±2.75) 


(-10.3, ±1.70) 


(-20.0, ±4.88) 


32 


(2.25,0.00) 


(-11.7, ±2.51) 


(-18.1,0.00) 


(-10.0, ±2.50) 


(-17.5,0.00) 


34 


(2.10,0.00) 


(-7.60,0.00) 


(-14.1,0.00) 


(-9.35,0.00) 


(-10.5,0.00) 



Fig. 8. — Contour plots of the </?— components of the fundamental mode, 6° , and of the first 
dipolar overtones, b^, b^, for model 4, model 29, model 31, and model 34 (in columns from 
left to right). The first overtone of model 4 (first column) and &^ of model 29 and model 
31 are complex and only their real part is shown. The aspect ratio and the local Rossby 
number increase from model 4 to model 34. Each contour plot is normalized by its maximum 
absolute value and the gray-scale coding varies from -1 (white) to +1 black. The contour 
lines correspond to ±0.1, ±0.3, ±0.5, ±0.7, and ±0.9. 
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converge to zero, in other words, the eigenvalues of the overtones approach that of the leading 
eigenmode. This appears to result from the fact that the structure of the fundamental 
mode becomes more similar to the following overtones, too. Differences between the first 
eigenmodes are most visible in their (y9-components. Figure [S] shows contour plots of the 
fundamental mode and the first two dipolar eigenmodes for model 4, model 29, model 31, 
and model 34. The axisymmetric azimuthal field of dynamo models at low Rossby number 
is dominated by two flux bund les of inverse pola rity close to the equatorial plane near the 
outer shell boundary (see also lOlson et al.lll999l ). These flux patches are visible in the 
components of the fundamental modes, 6°, for all four models. However, as the aspect 
ratio and the Rossby number increase, the axisymmetric flux concentration at low latitudes 
becomes less pronounced. The fundamental mode for model 34 resembles in this respect its 
following overtones and it is probably this adjustment in the eigenmodes which causes the 
convergence of the eigenvalues for dynamo models at Roe ^0.1. 

A systematic change in the field topology of the mean toroidal field with increasing 
Rossby number can be caused either by a change in the mean flow or by changes in the 
mean-field coefficients. The first possibility can be ruled out for dipolar models in the low 
Rossby number regime. For the sequence of models displayed in Figs. [7] and [H], the mean 
flow does not change significantly. Its inf luence on the dynamo mechanism is rela tively weak 
and not of primary importance (see also lOlson et al.lll999l : ISchrinner et al.l 120071 ). However, 
the mean-field coefficients, which are mainly responsible for the generation of the azimuthal 
field, indeed seem to vary in a systematic manner. Figure M displays the dominant dynamo 
coefficients acting on the toroidal field for model 4, model 29, model 31 and model 34. These 
are — 7^ and —751 describing an advection of the mean azimuthal field in radial and latitudinal 
direction, and —arr and —aeg leading to the generation of toroidal from poloidal field. We 



normalized all dynamo coefficients by the maximum modulus of for each model. This 
third diagonal component of the a-tensor is of crucial importance for the generation of 
poloidal from toroidal field and remains almost unchanged for models at different Ro£. 

The 7-effect advects the mean azimuthal field towards the equatorial plane and the 
outer shell boundary at low latitudes, and in the opposite direction at higher latitudes 
and in deeper layers (see also ISchrinner et al.l 120071 ). It is directly related to the columnar 
convection in a spherical shell. A column of fluid elements transported in an upwelling 
towards the outer spherical boundary has to shorten. Because of mass conservation, this 
causes a converging flow towards the equatorial plane. Vice versa, the rotation al constraint 



leads to a diverging flow in downwellings in deeper layers and at higher latitudes (jOlson et al. 
19991 ). For models approaching Roi ~ 0.1, the rotational constraint is relaxed and hence 
the 7-effect is diminished (see also Fig. [T0|) . In particular, the advective velocity towards 



the equatorial plane visualized by the outer contours of —'jg in Fig. [9] is less prominent for 
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Fig. 9. — Dynamo coefficients for model 4, model 29, model 31 and model 34 (in columns). 
The mean azimuthal field is advected in radial and latitudinal direction by — 7^ and — 7e, and 
generated from poloidal field by — a^r and —aee- For each model, the dynamo coefficients 
were normalized by the maximum modulus of a^i^. White stands for negative and black for 
positive values. The contour lines are =b0.075, ±0.05, and ±0.025 for aee and as in Fig. [8] 
for all other coefficients. 



-20- 



model 34 than for model 4. This is consistent with the changes in the topology of 6° for 
these models. 

The significance of the 7-effect in dipolar dynamo models is demonstrated by a mean- 
field calculation, in which 7 was arbitrarily suppressed. Figure [11] shows the azimuthally 
and time-averaged magnetic field for model 34 obtained from DNS (first row) and the lead- 
ing eigenmode of a corresponding mean-field calculation based on all mean-field coefficients 
determined (second row). Both are in good agreement. If the dynamo-coefficients related 
to the 7-effect are set to zero in a numerical experiment, the result changes substantially. 
The eigenvalue spectrum is flat and there are two complex, growing eigenmodes of either 
symmetry, i.e. the model is no longer kinematically stable. The resulting first real, dipolar 
eigenmode shown in the third row of Fig. [TT] varies on fairly small length scales. Moreover, 
the mean azimuthal flux concentrations at low latitudes near the outer boundary character- 
istic for dipolar dynamo models disappeared. 

Similar to the 7-components, a^r decreases considerably with increasing Rossby number, 
as demonstrated clearly by Fig. [TUJ It is the dominating coefficient responsible for the 
toroidal field generation by an a-effect. At low Rossby numbers, a strong 7-effect leads to 
the distinctive azimuthal field configuration with two flux bundles of inverse polarity close 
to the equatorial plane. The also increased Orr-component in this dynamo regime sustains it 
against the efficient diffusion due to strong gradients necessarily related to this field topology. 
The a^e-component, on the other hand, is on average much smaller than and remains at 
low level independent of the Rossby number. 
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Ro,=0.015 
Ro,=0.047 
Ro,=0.054 
Ro,=0.069 
Ro,=0.090 




Fig. 10. — The time and spatially averaged values of a^r, a^g, 7^, and 79 normalized by the 
mean value of for model 4, model 29, model 31, model 32, and model 34. 
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Fig. 11. — First row: Azimuthally and time-averaged magnetic field for model 34 obtained 
by direct numerical simulations. Second row: Leading eigenmode derived from ([7]) with the 
complete dynamo operator for model 34 as defined in ([8]). Third row: Eigenmode derived 
from ([7]) with a dynamo operator for model 34 in which the 7-effect was artificially sup- 
pressed. Each component is normalized by its maximum modulus and the gray-scale coding 
varies from -1 (white) to -|-1 (black). The contour lines correspond to ±0.1, ±0.3, ±0.5, ±0.7, 
and ±0.9. 
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6.2. Reason for the bistability and dynamo waves 



The existence of a dipolar and a multipolar dynamo regime with a sharp transition 
between both at Roe ~ 0.1 does not depend on the choice of the mechanical boundary 
conditions. The transition her e is therefore not controlled by the thickness of the Ekman 
boundary layers (as opposed to lKing et al.ll2009l ). However, most of the models with stress- 
free conditions at both boundaries and Roe < 0.1 exhibit a second, non-dipolar magnetic- field 
branch. Stress-free mechanical boundary conditions allow for the development of a strong, 
axisymmetric azimuthal flow, if the magnetic field is initially weak. In the presence of 
stress-free boundary conditions, the zonal geostrophic flow can only saturate owing to bulk 
viscosity and magnetic forces. When no-slip boundaries are present, the viscous braking of 
the geostrophic flow occurs mainly in boundary layers, an effe ct much stronger (by a factor 



than bulk viscous effects (e.g. iMorin &: Dormyl |2006| ). As a result, for stress-free 



conditions at the inner and the outer shell boundaries, even a very weak inertial forcing 
(Reynolds stresses) can yield a very significant geostrophic flow if the magnetic field is weak. 
The zonal flow pattern is then highly geostrophic, i.e. is constant on cylinders parallel 
to the rotation axis. Conversely, for the dipolar branch, the development of a zonal flow is 
partly inhibited by Maxwell stresses, le ading to substantial deviations from geostrophy (as 
also observed by lBusse fc Simitevll2009l ). The difference between both flow profiles is visible 
in Fig. [12] for model 45m and model 45d. Furthermore, the variation in the mean zonal flow 
causes differences in the f2-effect, rBr d{r~^V(j))/dr+r~^ sin 9 Bg d{sm 9~^V^)/d9, and thus 
in the dynamo mechanism for the toroidal field of both branches. For model 45m, the f2-effect 
correlates nicely with the mean azimuthal field (upper panel of Fig. [T2l) . and therefore, the 
model may be characterized as an afi-dynamo. However, for model 45d, the mean azimuthal 
field is for the most part not the result of the f2-effect. In particular the flux portions at 
larger radii and close to the equatorial plane are, if at all, anticorrelated to it. 

The difference in the dynamo mechanism for both branches explains why the Rossby 
number criterion for the dipolarity of the magnetic field applies only to the dipolar branch. 
The Rossby number criterion as formulated in this work compares the convective length scale 
with the Rossby radius. If the magnetic field is not solely a result of columnar convection but 
its generation mechanism also involves a large-scale zonal flow, the Rossby number criterion 
becomes meaningless. Other exaniples fo r whi ch it does not apply , presu mably for the same 
reason, can be found in lHori et al.l ( 2010 ) and Landeau fc Aubert ( 2011 ). Similar to models 
at high Ro£, model 45m lacks the particular azimuthal field configuration typical for dipolar 
or 'single-mode' dynamos. The magnetic field is governed on this branch by several modes 
and the relative dipole field strength drops to /dip = 0.25. The model is kinematically 
unstable and the magnetic field exhibits quasi-periodic time variations (see Fig. [5]). 
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Dynamo models of this stress-free multipolar branch and those with Roi > 0.1 are, for 
different reasons, not dominated by a dipolar mode. The dominance of only one, real dipolar 
eigenmode associated with the columnar flow is broken and in general complex overtones play 
an essential role in the dynamics of the magnetic field. If the magnetic Reynolds number is 
sufficiently low, the magnetic field evolves in the form of coherent dynamo waves. However, 
with increasing distance of the models from the dynamo threshold, the temporal coherence 
is lost presumably due to the enlarged number of relevant modes. A particular dynamo 
mechanism, on the other hand, is not the primary reason for oscillatory dynamos. Models 
35-3 8 (Roi > 0.1) exhibi t fairly coherent dynamo waves, but they are not of a aQ-tjpe 



see 



Schrinner et al.ll201ll ). For these examples, oscillatory dynamos are found because the 
rotational constraint is relaxed, i.e. Roe > 0.1, and nevertheless Rm remains moderate 
{Rm < 200). This twofold condition can be fulfilled for models with a thin convection 
zone, for exa mple, and governs the tran sition from steady to oscillatory dynamos already 
highlighted by iGoudard fc Dormvl ( 120081 ) . 



Christensen fc AubertI (120061 ) pointed out that the cube of the magnetic field strength 



for models with rigid boundaries is proportional to the measured buoyancy flux. The Lorentz 
number and the fiux-based Rayleigh number. Rag = Ra (Nu — l) E'^/Pr, can then be related 
trough 



Lo 

^1/2 
J ohm 



oc -Ra, 



1/3 



(9) 



where /ohm is the ratio of ohmic to total dissipation. The coefficient of proportionality was 
determined, still in th e case of rigid bou ndaries, to be 0.79 for dipolar models and 0.48 
for multipolar models ( jChristensenI 120101 ). We found, see Figil3l that the magnetic field 
strength for our limited sample of stress-free models is consistent with this relation, with 
similar prefactors. 



As noted already by iBusse fc SimitevI (120111 ) , there is almost no difference in the Nus- 
selt number, and thus in the flux-based Rayleigh number, for pairs of dipolar and multipolar 
models in the bistable regime (see Table [3]). However, /ohm is always smaller for the mul- 
tipolar branch, i.e. the afi-mechanism is less efficient than the a^-mechanism related to 
columnar convection. This deficiency together with the somewhat lower prefactor relates to 
for the lower field strength found for these models. 



6.3. Bifurcations between dynamo branches 



It is interesting to ponder on the transitions between the dipolar and multipolar branch 
for stress-free models when one single control parameter is varied. The two branches are 
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Fig. 12. — The axisymmetric azimuthal flow, the fi-effect expressed as rBrd{r~^V^)/dr + 
sin 6 Bg d{sm 9~^V ^) / dO , and the axisymmetric azimuthal magnetic fleld for model 45m 
(upper panel) and model 45d (lower panel). The style of the contour plots is explained in 
the caption of Fig. [TTl 



Table 3. Nusselt number and ratio of ohmic to total dissipation, /ohm, for bistable models 

with stress-free boundary conditions 











Nu and /ohm 


for models with stress-free boundary conditions 










Model 
Nu 

/ohm 


43m 
2.0 
0.29 


43d 
1.9 

0.48 


45m 
2.2 
0.28 


45d 46m 
2.4 3.0 
0.43 0.32 


46d 47m 47d 48m 48d 49m 49d 
3.0 3.0 3.1 2.9 2.9 2.3 2.5 
0.56 0.30 0.56 0.31 0.42 0.25 0.47 


51m 
1.3 
0.48 


51d 
1.4 
0.58 


52m 
4.3 
0.38 


52d 
4.2 
0.59 
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Fig. 13. — Lorentz number corrected by the dissipation factor versus the flux based Rayleigh 
number. Th e straight hnes rep resent the scahng for dipolar and multipolar models with rigid 
boundaries ( IChristensenI 120101 ). Our models with stress-free boundaries (squares stand for 
dipolar and diamonds for multipolar models) follow a similar scaling. 



illustrated in Fig. [14] for a fixed Ekman number of i^^ = 10~^ and magnetic Prandtl number 
of Pm = 1. For both branches, the local Rossby number increases with increasing Rayleigh 
numbers (see Table 1). If the Rayleigh number is increased from Ra = 1110 (model 52d) on 
the dipolar branch to Ra = 1200, the relative dipole field strength collapses (the local Rossby 
number crosses the Roi ~ 0.1 boundary). The multipolar field configuration then appears 
to be the only stable solution (circle on the figure) and a hysteretic behavior is observed if 
the Rayleigh number is decreased from this state (i.e. model 52m is then obtained). 

Interestingly, the transition between both branches is not always as abrupt as in Fig. 
[m Instead, the two branches can also merge more continuously if the zonal geostrophic 
flow on the multipolar branch becomes too weak. This is best demonstrated by varying the 
magnetic Prandtl number keeping all other control parameters fixed. Figure [T5] presents 
the axisymmetric toroidal kinetic energy as a function of the magnetic Prandtl number at 
fixed Rayleigh number {Ra = 750) and Ekman number {E = 10~^). Both branches are 
represented (the multipolar solution corresponding to the larger values). If the magnetic 
Prandtl number is increased the saturated value of the geostrophic flow on the multipolar 
branch decreases. For Pm = 6 ( circle on the figure) the multipolar solution is only meta- 



stable (see iMorin &: Dormyl |2009| ) . It can be observed for a short period of time (enough to 
assess its amplitude) but then switches to the dipolar solution. For this value of the control 
parameter (Pm), only the dipolar solution could be produced. The saturated amplitude of 
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the geostrophic flow is here too small to prevent the growth of the dipolar solution. This 
behavior is typical for simulations at large Pm, although the value at which the multipolar 
branch is lost depends on the other parameters, too. 

Finally, besides the transition between both branches, it is worth pondering on how 
these dynamo solutions bifurcate from the purely hydrodynamic solution. We have not 
performed a detailed study of this problem, but we observed in Fig. [15] that for Pm = 
0.5, both solutions exist, while the hydrodynamic solution is stable (i.e. a small magnetic 
perturbation decays, and a finite amplitude perturbation is needed to obtain any of the 
dynamo sol utions). The subc r itical bifurcation of the dipolar mode at low Pm was already 



described in iMorin fc Dormyl (120091 ) with no-slip boundary conditions. It is interesting that 
the system exhibits here a triple stability (two dynamo solutions and a purely hydrodynamic 
mode). Obviously the dynamo bifurcation with stress-free boundary conditions deserves 
further studies. 



-28- 



<1.5- 



@ 



1000 1500 

Ra 



Fig. 14. — Evolution of the magnetic field strength, measured by the Elsasser number, 
for both branches as the Rayleigh number is varied at fixed Ekman and magnetic Prandtl 
numbers {E — 10~^ and Pm — 1). 




Fig. 15. — Axisymmetric toroidal energy density of the flow as a function of the magnetic 
Prandtl number at fixed Rayleigh number {Ra = 750) and Ekman number {E = 10^^) 
defined as one half of the average over the fluid volume of the non-dimensional toroidal 
velocity squared. 
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6.4. Dynamo models versus stellar and planetary dynamos 



It is not known to which extent simplified dynamo models indeed refiect physical pro- 
cesses going on in stellar or planetary dynamos. All simulations to date are performed in 
a wrong parameter regime; they are therefore not directly comparable with observations. 
However, scaling laws derived from numerical dynamo simulations could help to test their 
reliability. Results obtained from dynamo models may be compared with the strength, the 
geometry and the time dependence of stellar and planetary magnetic fields. 

The relation between the fi e ld str ength and the fiux-based Rayleigh number (Q pro- 
posed by IChristensen fc AubertI (120061 ) is co nsistent with the field-st rength for a class of 
fast rotating stars and some of the planets (IChristensen et al.l l2009l ). However, without 
further assumptions, it is neither apphcabl e to slowly rotati ng stars, e.g. the Sun, nor 
to Mercury, Saturn, Uranus and Neptune (jChristensenI 120101). For theses e xamples, the 
field strength falls below the predicted value. Moreover, iMorin et al.l (|2010( ) observed M 
dwarfs of very similar mass (i.e. energy fiux) and rotation rate, but with dyna mo generated 



magne tic fields which differ in their field strength and their field topologies. iMorin et al. 



( I2OIII ) argued that a different force balance could be responsible for the observed bistabil- 
i ty, similar to the strong-field branch scenario introd uced in the context of the geodynamo 
( iRoberts. P. H. fc Soward. A. M.lll978l : lRobertslll988l ). In this study, we point out that differ- 
ences in the dynamo mechanism could likewise lead to two different magnetic field branches. 
The shearing of poloidal field lines (fi-effect) due to a strong geostrophic zonal fiow may 
play an important role for the field generation in multipolar, but not in dipolar models. A 
dynamo mechanism based on the action of a mean zonal fiow leads to models characterized 
by a more variable and a somewhat weaker magnetic field. The bistability investigated in 
this study was found for a wide range of parameter values. It could also account for different 
magnetic fields of stars with similar parameters at a Rossby number close to 0.1, which has 
not been observed yet. 

The Rossby number criterion for the dipolarity of the magnetic field is co nsistent with 

the t opology of planetary magnetic fields, except for Uranus and Neptune fsee IChristensen 

2010). Furthermore, it gets some support from observations of stellar magnetic fields (IMorin et al 
20081 ). In particular, a decrease of the size of the convection zone may lead to higher local 
Rossby numbers and th us to an abrupt change in the field topology, as observed for early M 



dwarfs (see Fig. l2land lMorin et al.ll2008l : iDonati et al.ll2008 



Reiners fc Basril 120091 ) . Like- 



wise, the multipolar magnetic fields of Uranus and Neptune could be compatible with the 
Rossby numb er rule, if convection in these pla nets takes place in a thin convection zone, as 
suggested by Istanlev fc BloxhanJ J2004j . bood ). 



The regime boundary at Roi ~ 0.1 separates models with a stable dipole field from 
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models exhibiting dipole reversals. Why are models of the dipolar regime non-reversing? 
The fact that modes, other than the fundamental mode, are on average damped does not 
prevent reversa ls, which may be tr iggered by the coupling with a weakly damped competing 
dynamo mode (jPetrelis et al.ll2009l ). But the 'single mode property' of these models reported 
above excludes that different modes become critical and explains why a reversal mechanism 
based on the coupling of competing modes has not been identified in these models. It has 
been proposed that the value of Ro^ for the Earth's core may be about 0.1 and thus in the 
vicinity of the regime bound ary, which could also provide a mechanism for polarity reversals 
doison fc Christensen|[2006[ ). 



Whether stellar and planetary magnetic fields are dipolar for the same reason as the 
models could perhaps be assessed studying their time dependence. The single mode property 
should lead to differen t time scales for the variation of the dipolar and the nondipolar field 



(ISchrinner et al. 



Le Mouellll984[ ). 



20111). as it has been repor ted for the geomagnetic secular variation (e.g. 



Tanriverdi fc Tilgnerl ( 1201 ll ) recently pointed out that single mode dynamos 



may be identified by the spectrum of temporal fluctuations of the magnetic energy. If the 
velocity spectrum is characterized by white noise and the evolution of the magnetic field 
is dominated by only one single dynamo mode, a o;"^ dependence for the spectrum of the 
magnetic field was predicted. This might be verifiable with magnetic field data. As discussed 
above, reversals of the Earth's magnetic field clearly indicate that the geodynamo cannot 
always match the single mode property. 

We found coherent magnetic cycles in our non- dipolar dynarno mod els only if the mag- 
netic Reynolds number is sufficiently low. Similarly, iBrown et al.l ( 1201 ll ) reports oscillations 
of the magnetic field in an anelastic simulation at fairly low magnetic Reynolds number, 
Rm = 136. If Rm is increased in our simulations, the temporal coherence is lost. Moreover, 
the dynamo waves in our models migrate from low latitudes towards the poles (see also 
Schrinner et al.l 1201 if ), as opposed to solar sunspot regions. Given that estimates for the 
magnetic Reynol ds number range f rom 10^ in the photosphere to 10^° at the base of the 
convection zone (lOssendrij verl |2003| ) . it remains unclear how the temporal cohe rence visible 



in th e 22-year solar cycle can persist in such a highly turbulent environment ( iJones et al. 
2010h . 



7. Conclusions 

Convection in a rapidly rotating spherical shell is organized in quasi-geostrophic columns 
parallel to the rotation axis. It gives rise to highly efficient but also very selective dynamo 
action: only one, real, dipolar eigenmode of the magnetic field is sustained. This single-mode 
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property accounts for the dipole dominance and the stabihty of the dipole field in models 
of this regime. Consequently, the dipolarity of our models collapses, if the dominance of 
the fundamental mode is broken and, in general, complex overtones dominate the evolution 
of the magnetic field. This may happen if convection is less constrained by rapid rotation, 
or, if the magnetic field is not solely a result of columnar convection, but a mean zonal 
flow strongly influences the dynamo mechanism. Whether or not differences in the topology 
and the time variability of planetary and stellar magnetic fields may be explained by the 
dichotomy between these two dynamo regimes certainly needs to be further explored. 

We acknowledge discussions with J.F. Donati and J. Morin. MS is grateful for finan- 
cial support from the DFG fellowship SCHR 1299/1-1. LP and ED acknowledge financial 
support from "Programme National de Physique Stellaire" (PNPS) of CNRS/INSU, France. 
Computations were performed at CINES and CEMAG computing centers. 



A. Appendix 

The definition of the local Rossby number given here relies on a length scale de- 
rived from the mean harmonic degree of only the poloidal velocity field. This intends 
to filter out the contribution of the mean toroidal flow, which is not negligible in mod- 
els with_jtress;ft;ee_jn^clm conditions. The local Rossby number introduced 



by IChristensen fc AubertI (120061 ) based on a length scale derived from the mean harmonic 
degree of the total flow does not allow us to distinguish dipolar models from multipolar 
models equally well (Fig. Uni left panel). 

In this study we argued that the typical length scale of convection relative to the Rossby 
radius strongly influences the topology of the magnetic field in our models. The definition 
of a local Rossby number based on a velocity field, in which only the axisymmetric toroidal 
contribution is canceled out is maybe more appropriate to test this argument than the 
definition used throughout in the text. The right panel of Fig. [16] shows that a local Rossby 
number defined in this way seems to serve equally well to distinguish dipolar from multipolar 
models. 
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Fig. 16. — Relative dipole field strength versus local Rossby number. On the left, the local 
Rossby number is base d on a length scale derived fro m the mean harmonic degree of the total 
flow as introduced by IChristensen &: AubertI (120061 ) . on the right, it is based on a velocity 
field and its typical length scale without considering the contribution of the mean toroidal 
flow component. The meaning of the symbols is the same as in Fig. [H 
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